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Abstract. In our model with a Standard Model gauge group extended with a baryon 
number minus lepton number charge for each family of quarks and leptons, we calculate 
the baryon number relative to entropy produced in early Big Bang by the Fukugita- 
Yanagida mechanism. With the parameters, i.e., the Higgs VEVs already fitted in a 
very successful way to quark and lepton masses and mixing angles we obtain the order 
of magnitude pure prediction Yb = 2.59 ^j^j^ 10~^^ which according to a theoretical 
estimate should mean in this case an uncertainty of the order of a factor 7 up or down 
(to be compared to Yb = (1.7 — 8.1) x lO"^'^) using a relatively crude approximation for 
the dilution factor, while using another estimate based on Buchmiiller and Pliimacher 
a factor 500 less, but this should rather be considered a lower limit. With a realistic 
uncertainty due to wash-out of a factor 100 up or down we even with the low estimate 
only deviate by 1.5ct. 



1 Introduction 

Using the model for mass matrices presented by us in an other contribution [Q| 
at this conference we want to compute the amount of baryons produced in the 
early universe. This model works by having the mass matrix elements being 
suppressed by approximately conserved quantum numbers from a gauge group 
repeated for each family of quarks and leptons and also having a {B — L) charge 
for each family. 

The baryon number density relative to entropy density, Yb, is one of the 
rather few quantities that can give us information about the laws of nature 
beyond the Standard Model and luckily we have from the understanding of 
the production of light isotopes at the minute scale in Big Bang fits to this 
quantity |^] . The "experimental" data of the ratio of baryon number density to 
the entropy density is 



= (1.7-8.1) X 10"" . (1) 

cxp 



We already had a good fit of all the masses and mixings [^|JJ] for both quarks and 
leptons measured so far and agreeing with all the bounds such as ncutrinolcss 
beta decay and proton decay not being seen and matching on the borderline but 
consistent with the accuracy of our model and of the experiment of CHOOZ the 
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electron to heaviest left-handed neutrino mixing, and that in a version of our 
model in which the dominant matrix element in the right-handed neutrino mass 
matrix is the diagonal one for the "third" {i.e. with same (B — L)i as the third 
family) family j/i{3 right-handed neutrino. This version of our model which fits 
otherwise very well does not give sufficient [B — L) excess, that survives, but the 
by now the best model in our series should have the right-handed mass matrix 
dominated by the off-diagonal elements (2,3) and (3,2), so that there appears 
two almost mass degenerate see-saw neutrinos, in addition to the third one (first 
family) which is much lighter. 



2 Mass matrices and results for masses and mixing angles 

Our model produces mass matrix elements - or effective Yukawa couplings - 
which are suppressed from being of the order of the top-mass because they are 
forbidden by the conservation of the gauge charges of our model and can only 
become different from zero using the 6 Higgs fields [Q,^ which we have in addition 
to the field replacing the Weinberg-Salam one. In the neutrino sector according 
to the see-saw mechanism |^ we have to calculate Dirac- and Majorana-mass 
matrices, Meg ~Af^M^^ (M^)^, to obtain the effective mass matrix M^e for 
the left handed neutrinos we in practice can "see" . Here we present all mass 
matrices as they follow from our choice of quantum numbers for the 7 Higgs 
fields in our model and for the quarks and leptons (as they can be found in the 
other contribution). Only the quantum numbers for the field called (pB-L is - 
in order to get degenerate see-saw neutrinos - changed into having the B — L 
quantum numbers of family 2 and 3 equal to 1, i.e., {B — L)2 = (B — L)^ = 1, 
while the other family quantum numbers are just zero: 
the up-type quarks: 



M„ 



((0ws)t) 



V2 



the down-type quarks 
the charged leptons: 



V2 



(wt)3iytr2 ^ptv^tT2 ^p^w'f fT\ 



(2) 



(3) 



w3iy(Tt)2 (uj^fp^W{T^f {uo^fp^WT^^ 
LO^{p^fW(T^f W{T^f WT^X 

V2 \w6(pt)3(^t)2r4 (^1)2^4 



the Dirac neutrinos 



, (pt)3w^tT2 iytT2 
V2 \_(pt)3yytTt^t lytyt^t 



(4) 



(5) 
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and the Majorana (right-handed) neutrinos: 

/ (pt)6^t (pt)3^t/2 (pt)3/2\ 

Mfl,^(03-.> (P^)'xV2 1 (6) 

V(P^)V2 1 X I 

We shall remember that it is here understood that all the matrix elements are 
to be provided with order of unity factors which we do not know and in practice 
have treated by inserting random order of unity factors over which we then 
average at the end (in a logarithmic way). 



3 Renormalisation group equations 

The model for the Yukawa couplings we use gives, in principle, these couplings 
at the fundamental scale, taken to be the Planck scale, at first, and we then use 
the renormalisation group to run them down to the scales where they are to be 
confronted with experiment. From the Planck scale down to the see-saw scale 
or rather from where our gauge group is broken down to SAIG x U{1)b-l we 
use the one-loop renormalisation group running of the Yukawa coupling constant 
matrices and the gauge couplings 1^ in GUT notation including the running of 
Dirac neutrino Yukawa coupling: 

.(,2^91 41 3 2^92 19 3 ,(.,2^93 » 3 



_ ,dY^ 3 
Ibn —— = - 
dt 2 



(Y^iY^)^ -YuiYu)^) + Q.g? + I9I + 833') I Y^, 

{Y4Y.)^ - Y{Y)^) Y. + i^Ys- (^g? + I9I) } 



1671^^ = ^ (YJY^V - YAY^V) Y, + iYs-[ -q1 + -fl,^ I \ Y, 

_ ^dY^ 3 
dt 2 



(K(K)t _ y^(y^)t) + I - (Ag? + \gi^ | K , 

Ys =Tr(3yJy„ + 3y^^y^ + rir^ + ytK) , 

where t = In/z and /i is the renormalisation point. 

In order to run the renormalisation group equations down to 1 GeV, we use 
the following initial values: 



U{1): 


9i{Mz) = 


0.462 , 


5l(Afpianck) 


= 0.614, 


(7) 


SU{2) : 


92{Mz) = 


0.651 , 


52(Afrlanck) 


= 0.504, 


(8) 


SU{3) : 


93{Mz) = 


= 1.22 , 


53(-A^Planck) 


= 0.491 . 


(9) 



We varied the 6 free parameters and found the best fit, corresponding to the low- 
est value for the quantity g.o.f. = ^ In f -^7^7^^) = 3.38, with the following 
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values for the VEVs: 

{(j)„s) = 246 GeV , {(Pb-l) = 1-23 x 10^" GeV , {u) = 0.245 , 
{p) = 0.256 , {W) = 0.143 , (T) = 0.0742 , (x) = 0.0408 , (10) 

where, except for the Weinberg-Salam Higgs field and {(j)B-L), the VEVs are 
expressed in Planck units. Hereby we have considered that the Weinberg-Salam 
Higgs field VEV is already fixed by the Fermi constant. The results of the best 
fit, with the VEVs in eq. (|l^), are shown in Table |l|. 

4 Quantities to use for baryogenesis calculation 

Since the baryogenesis in the Fukugita-Yanagida scheme arises from a nega- 
tive excess of lepton number being converted by Sphalerons to a positive baryon 
number excess partly and this negative excess comes from the CP violating 
decay of the see-saw neutrinos we shall introduce the parameters giving the 
measure of the relative asymmetry under C or CP in the decay of neutrino 
number i: Defining the measure for the CP violation 

where P are N^i decay rates (in the Nui rest frame), summed over the neutral 
and charged leptons (and Weinberg-Salam Higgs fields) which appear as final 
states in the Nj^i decays one sees that the excess of leptons over anti-leptons 
produced in the decay of one N^i is just e^. The total decay rate at the tree level 
is given by 

Pn. = Pn., + PNa = ^V/i ^2 > (12) 

47r {(pivs) 





Fitted 


Experimental 




Fitted 


Experimental 


mu 


5.2 MeV 


4 MeV 


rrid 


5.0 MeV 


9 MeV 


nic 


0.70 GeV 


1.4 GeV 


rus 


340 MeV 


200 MeV 


Mt 


208 GeV 


180 GeV 


rub 


7.4 GeV 


6.3 GeV 


me 


1.1 MeV 


0.5 MeV 




0.10 


0.22 




81 MeV 


105 MeV 




0.024 


0.041 


rrir 


1.11 GeV 


1.78 GeV 




0.0025 


0.0035 


ArriQ 


9.0 X 10"'' eV^ 


4.5 X lO"'' eV^ 




1.8 X 10-" eV^ 


3.0 X lO-" eV^ 


tan^ 9q 


0.23 


0.35 


tan &atm 


0.83 


1.0 


tan ^chooz 


3.3 X 10"^ 


< 2.6 X 10"^ 


g.o.f. 


3.38 





Table 1. Best fit to conventional experimental data. All masses are running masses 
at 1 GeV except the top quark mass which is the pole mass. Note that we use the 
square roots of the neutrino data in this Table, as the fitted neutrino mass and mixing 
parameters (m), in our goodness of fit (g.o.f.) definition. 
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where can be expressed through the unitary matrix diagonalising the right- 
handed neutrino mass matrix Vf>: 



MP 



Mr Ml Vr = diag ( Ml , M| , Ml 
The CP violation rateis computed according to 



Y.,^^iHmPVMp)%] f 



(13) 
(14) 



(15) 



where the function, f{x), comes from the one-loop vertex contribution and the 
other function, g{x), comes from the self-energy contribution. These e's can be 
calculated in perturbation theory only for differences between Majorana neutrino 
masses which are sufficiently large compare to their decay widths, i.e., the mass 
splittings satisfy the condition, \Mi — Mj\ ^ {Fi — Fjl: 



1 - (1 + x) In 



1 + x 



l-x 



(16) 



We as usual |g] introduce the dacay rate relative to 

Mpianck {{MpyMP), 



K, 



2H 



T=Mi 



1.66 < 



M,, 



(z = 1,2,3) , (17) 



where Fi is the width of the flavour i Majorana neutrino, is its mass and 
g^i is the number of degrees of freedom at the temperature Mi (in our model 
~ 100). 

In order to estimate the effective K factors we first introduce some normalized 
state vectors for the decay products: 



E| [M,?(Af.) 



ki 



MP [Mi) _ , MP{M, 



MP{M,) 



3i 



Then we may take an approximation for the effective K factors: 
K.„i=^i(Mi) , 

ifo«2 = K2{M2) + |(2|3)|' i^3(M3) + |(2|1)|' Ki{Mi) , (20) 



i^3(M3) + |(3|2)|' K2{M2) + |(3|1)|' Ki{Mi) . 



(18) 

(19) 
(20) 
(21) 



5 Result for baryogenesis 



Using the Yukawa couplings - as coming from the VEVs of our seven different 
Higgs fields ~ the numerical calculation of baryogenesis were performed using 
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our random order unity factor method. In order to get baryogenesis in Fukugita- 
Yanagida scheme, we calculated the see-saw neutrino masses, i^ctti factors and 
CP violation parameters using N = 10, 000 random number combinations and 
logarithmic average method: 



Ml 


^ 2.1 X 10^ 


GeV 


, K^tti 


= 31.6 , 


kil 


= 4.62 X 10 


-12 


Ms 


= 8.8 X IQS 


GeV 




= 116.2 . 


.hi 


= 4.00 X 10 


-6 


Ms 


= 9.9 X IQS 


GeV 




= 114.7 . 


.hi 


= 3.27 X 10 


-6 



The sign of ei is unpredictable due to the complex random number coefficients 
in mass matrices, therefore we are not in the position to say the sign of e's. Using 
the complex order unity random numbers being given by a Gaussian distribution 
we get after logarithmic averaging using the dilution factors as presented by |^,|| 

ys = 2.59™ ^ 10"" ' (22) 

where we have estimated the uncertainty in the natural exponent according 
ref. to be 64 % • VTO « 200 %. 

The understanding of how this baryon to entropy prediction Yb comes about 
in the model may be seen from the following (analytical) estimate 

^3 ^1.10-9 (23) 

3 ^/gZT^ Mpianck 3 

where we left out for simplicity the In K factor in the denominator of the dilution 
factor K and where il/3 is the mass of one of the heavy right-handed neutrinos in 
our model M3 w Since the atmospheric mass square difference square 

root \J Aml^tm ~ 0-05 eV ~ {4'ws) [WT^ / we see that keeping it leaves us 
with the dependence 

3v/0.05 eV- g.Mpunc^W^T^ ^ 5 ^*W^T^ ^ ' 

6 Problem with wash-out effects? 

To make a better estimate of the wash-out effect we may make use of the calcu- 
lations by by putting effective values for the see-saw neutrino mass M and 
rh. The most important wash-out is due to "on-shell" formation of right-handed 
neutrinos and only depends on K or the thereto proportional m, but next there 
are wash-out effects going rather than by K or m as Mfh? . In the presentation 
of the results by fixed ratios between right-handed neutrino masses were 
assumed. However, in reality a very important wash-out comes form the off-shell 
inverse decay and that goes as 

MiV^™,^ with ^,.K^1!^ (25) 
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Here we use the notation with from \11 : rrij Kj ■ 2.2 • 10 "^eV. 

Using such a term (see eq. ^) with the ansatz ratios used in |11 , M| = 
10^ and M| = 10^ one gets for eq. (|5j) « 10^ Mi fhl, while we would 
with our mass ratios (eq. |2|) M| w 1/4- 10^'' Mf and Af| « 1/4- 10^° obtain 
correspondingly 2 • 10^ GeV • 1/4 • 10^" m§ « 1/2 • lO^^ GeV mj, which then 
being identified with 10^ Mi use fni would lead to that we should effectively 
use for simulating our model the mass of the right handed neutrino - which 
is a parameter in the presentation of the dilution effects in - Mi use = 
1/2 • 10^5 GeV/lO^ = 1/2 • 10^ GeV. Inserting this Mi use value for our estimate 
m2 ~ m3 Ri 0.1 eV gives a dilution factor k w 10~^, i.e., a factor 500 less than 
what we used with our estimate using the K^s's. (Our fh^ = m2 are surprisingly 
large compared to the ^/Arn^^ because of renormalzation running .) Using the 
better calculation of Q which has a very steep dependence - a fourth power 
say - as function of fh our uncertainty should also be corrected to a factor 100 
up or down. So then we have one and a half standard deviations of getting too 
little baryon number. 



7 Conclusion 

We calculated the baryon density relative to the entropy density - baryogenesis 
- from our model order of magnitudewise. This model already fits to quark and 
lepton masses and mixing angles using only six parameters, vacuum expectation 
values. We got a result for the baryon number predicting about a factor only 
three less than the fitting to microwave background fluctuations obtained by 
BuchmuUcr et al. p^ , when we used our crude ifcff's approximation. However, 
using the estimate extracted from the calculations of |jll| we got three orders of 
magnitude too low prediction of the baryon number. This estimate must though 
be considered a possibly too low estimate because there is one scattering effect 
that is strongly suppressed with our masses but which were included in that 
calculation. But even the latter estimate should because of the steep dependence 
of the result on the parameters be considered more uncertain and considering 
the deviation of our prediction only 1.56(t is not unreasonable. 

Since we used the Fukugita-Yanagida mechanism of obtaining first a lepton 
number excess being converted (successively by Sphalerons) into the baryon 
number, our success in this prediction should be considered not only a victory 
for our model for mass matrices but also for this mechanism. Since our model 
would be hard to combine with supersymmetry - it would loose much of its 
predictive power by having to double the Higgs fields - we should consider it in 
a non SUSY scenario and thus we can without problems take the energy scale 
to inflation/reheating to be so high that the plasma had already had time to 
go roughly to thermal equilibrium before the right-handed neutrinos go out-of- 
equilibrium due to their masses. We namely simply have no problem with getting 
too many gravitinos because gravitinos do not exist at all in our scheme. 

Another "unusual" feature of our model is that the dominant contribution to 
the baryogenesis comes from the heavier right-handed neutrinos. In our model it 
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could be arranged without any troubles that the two heaviest right-handed neu- 
trinos have masses only deviating by 10% namely given by our VEV parameters 
X- This leads to significant enhancement of the 62 and 63 which is crucial for the 
success of our prediction. There is namely a significant wash-our taking place, 
by a factor of the order of k = 10~^ to 10~^. It is remarkable that we have here 
worked with a model that order of magnitudewise has with only six adjustable 
parameters been able to fit all the masses and mixings angles for quarks and 
leptons measured so far, including the Jarlskog CP violation area and most im- 
portantly and interestingly the baryogenesis in the early Universe. To confirm 
further our model we are in strong need for further data - which is not already 
predicted by the Standard Model, or we would have to improve it to give in prin- 
ciple accurate results rather than only orders of magnitudes. The latter would, 
however, be against the hall mark of our model, which precisely makes use of 
that we can guess that the huge amount of unknown coupling constants in our 
scheme with lots of particles can be counted as being of order unity. 
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